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The  purpose  of  this  analytical  study  was  to  develop  new  or  improved  altitude  ignition  and 
lean  deceleration  blowout  correlations  which  relate  the  operating  limits  of  a  combustion  system 
to  engine  operating  parameters.  These  correlations  can  then  be  universally  applied  to  current  or 
future  aircraft  gas  turbine  combustion  systems  for  ignition  and  stability  assessment. 

The  data  base  used  in  constructing  the  models  consisted  of  the  substantial  body  of 
experimental  data  acquired  during  recent  USAF,  Navy,  NASA  and  P&W  sponsored  programs  on 
operability  performance  of  military  and  commercial,  gas  turbine  engines  including  F100(3), 
PW1128  Development,  F100-ILC  Core,  PW2037,  JT9D-7R4  and  TF30-P-414/414A. 

The  approach  taken  in  this  study  was  to  consider  the  effects  of  the  following  physical 
processes  in  the  combustor  fundamental  to  ignition  and  flame  stabilization:  thermal  and 
turbulent  diffusion,  chemical  reaction,  fuel  evaporation  and  turbulent  mixing.  Since  more  than 
two  physical  processes  are  involved,  the  models  for  ignition  and  lean  deceleration  blowout  were 
formulated  based  on  the  characteristic  “time”  rather  than  “rate”  approach.  The  characteristic 
time  model  does  not  attempt  to  analyze  the  entire  combustor  flow  field,  but  rather  considers  key 
regions  of  the  flow  field  and  defines  expressions  for  the  time  required  for  spark  kernel  quenching, 
chemical  reaction,  fuel  evaporation  and  hot  flow  residence. 

The  criteria  for  successful  ignition  and  flame  stabilization  were  defined  as  follows: 

•  The  limiting  case  for  ignition  occurred  when  the  time  for  quenching  the 
spark  kernel  is  on  the  order  of  the  sum  of  the  chemical  reaction  plus 
evaporation  times  and  evaluated  in  the  region  near  the  igniter 

•  The  flame  stabilization  was  viewed  as  the  competition  between  the  hot  flow 
residence  time  in  the  primary  zone  and  the  chemical  reaction  plus 
evaporation  times  evaluated  in  the  shear-layer  region  between  the  hot 
recirculation  zone  and  the  free  stream. 

Table  1  presents  the  characteristic  times  which  have  been  related  and  the  combustor 
physical  processes  which  were  found  to  provide  the  best  models  for  altitude  ignition  and  lean 
deceleration  blowout.  The  first  column  in  Table  1  lists  the  performance  parameter  and  the  mode 
of  engine  operation  which  were  evaluated.  The  second  column  lists  the  combustor  limits  which 
were  determined  in  this  study.  The  third  column  lists  the  combustor  related  parameters  which 
were  determined  according  to  the  geometric  and  operating  characteristics  of  the  individual 
combustors.  The  fourth  column  lists  the  governing  physical  process  for  the  specific  models  which 
were  developed  in  this  study  and  the  last  column  identifies  the  degree  of  the  relationship  based 
on  the  regression  coefficient  of  determination. 

The  combustor  physical  processes  (evaporation  and  chemical  reaction)  and  operating 
parameter  (spark  kernel  quench  time)  developed  by  Ballal  and  Lefebvre,  Reference  1  for  the 
phenomena  of  spark  ignition  were  selected  to  predict  altitude  ignition  performance.  Since 
ignition  and  flame  stabilization  are  similar  phenomena,  the  combustor  physical  processes  for  the 
phenomena  of  flame  stabilization  were  defined  in  the  same  manner  as  in  the  spark  ignition  with 
the  exception  that  the  parameters  involved  are  evaluated  according  to  the  primary  zone 
conditions.  These  physical  processes  and  the  combustor  primary  zone  volumetric  hot  residence 
time  were  selected  to  predict  altitude  stability  and  lean  deceleration  blowout  limit. 


TABLE  1. 


ALTITUDE  IGNITION  AND  LEAN  DECELERATIONS  MODELS 


Performance 

Parameter 

Combustor 

Limit 

Combustor 
Operating  Parameter 

Governing  Physical 
Process 

Model 

Rating 

1.  Altitude  Ignition 

a.  Windmill 

Ignition 

Stability 

Spark  Kernel  Quenching 
Volumetric  Hot  Residence* 

Fuel  Evaporation 
Chemical  Reaction  and 
Fuel  Evaporation 

Very  Good 
Good 

b.  Spooldown 

Ignition 

Stability 

Spark  Kernel  Quenching 
Volumetric  Hot  Residence 

Fuel  Evaporation 
Chemical  Reaction  and 
Fuel  Evaporation 

Good 

Very  Good 

2.  Lean  Deceleration 
Blowout 

Stability 

Volumetric  Hot  Residence 

Chemical  Reaction 

Excellent 

•Combustor  Primary  Zone  Hot  Residence  Time 

148SC 


Although  the  evaporation  time  was  found  to  be  longer  than  the  chemical  reaction  time  for 
the  phenomena  of  windmill  and  spooldown  ignition,  consideration  of  both  physical  processes 
provided  a  higher  degree  of  correlation.  However,  the  lean  deceleration  blowout  limit  is 
controlled  by  chemical  reaction  for  the  engine  data  used  in  this  study. 

In  most  cases,  very  good  data  correlations  with  near  zero  intercept  have  been  obtained. 
Examination  of  the  correlations  provides  good  insight  into  the  effect  of  combustor  and  engine 
parameters  on  altitude  ignition  and  lean  deceleration  blowout  performance.  The  accuracy  of  the 
ignition  and  flame  stabilization  models  is  limited  by  the  inadequate  knowledge  of  the  mean  drop 
size,  drop  size  distribution,  turbulence  intensity,  and  equivalence  ratio  in  the  combustion  zone 
and  near  the  ignitor.  Further  improvements  in  the  models  will  require  measurements  of  drop  size 
and  turbulence  intensity  levels  to  better  quantify  these  variables. 
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(30)  (33) 

Characteristic  time  for  fuel  evaporation  (stability),  msec  Equations  (20)  (24) 

(31)  (34) 

Characteristic  time  for  chemical  reaction  (stability),  msec  Equations  (20)  (32) 
(35) 

Equivalence  ratio 

Fraction  of  total  fuel  in  form  of  vapor 


Subscripts 
a  Air 

comb  Combustor 

f  Fuel 

g  Gas 

i  Ignition 

s  Stability 

inj  Injected 

ign  Plane  of  the  ignitor 

o  Reference  value 

ad  Adiabatic  value 

st  Stoichiometric  value 

Primary  zone 
Combustor  inlet  value 


pz 
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Conversion  Factors 


To  Convert  From: 

To: 

Multiply  by: 

Length 

Centimeter  (cm) 

Inches  (in.) 

3.937  E-01 

Meter  (m) 

Feet  (ft) 

3.281 

Micron  (pm) 

Meters  (m) 

1.000  E-06 

Mass 

Kilogram  (kg) 

Pounds  (Ibm) 

2.205 

Force 

Newton  (kg-m/s2) 

Pounds  (lbf) 

2.248  E-01 

Pressure 

Newton/square  meter 
(N/m2) 

Atmosphere  (atm) 

9.869  E-06 

Kilopascal  (kPa) 

Pounds/square  foot 
(lbf/ft2) 

2.089  E+01 

Temperature 

Degree  Kelvin  (K) 

Degree  Fahrenheit  (°F) 

tp  «  1.8  [t|( 

Degree  Rankine  (°R) 

Degree  Fahrenheit  (°F) 

tR  =  tF  +  4. 

Degree  Kelvin  (K) 

Degree  Rankine  (°R) 

“  Wl.8 

Volume 

Cubic  meter  (m:i) 

Cubic  feet  (ft1') 

3.531  E+01 

Density 

Kilogram/cubic  meter 
<  kg/nr1) 

Pounds- mass/cubic  feet 
Obm/ft3) 

6.243  E-02 

Surface  Tension 

Newton/meter  (N/m) 

Pounds-force/foot  (lbf/ft) 

6.852  E-02 

Dynamic  Viscosity 

Kilogram/meter-second 

(kg/m-s) 

Pounds-force/foot-second 

(Lbf/ft-s) 

6.721  E-01 

Kinematic  viscosity 

Squared  meter/second 
(m2/s) 

Square-feet/second  (ft2/s) 

1.076  E+01 

Thermal  conductivity 

Joule/meter -second-Kelvin 
(J/m-s/K) 

Btu/hour-foot-degree 

Rankine 

(Btu/hr-ft-°R) 

5.780  e-01 

Specific  Heat 

Joule/kilogram-Kelvin 

(J/kg-K) 

Btu/pounds-mass-degree 

Fahrenheit 

(Btu/lbm-<>F) 

2.388  E-04 

Energy 

Joule 

British  thermal  units 
(Btu) 

9.486  E-04 

Velocity 

Centimeter/second  (cm/s) 

Inches/second  (m/s) 

3.937  E-01 

Meter/second  (m/s) 

Feet/second  (ft/s) 

3.281 

Physical  Constants 

Universal  gas  constant  Ru  =  0.730  atm-ft3/(lb-mol)  (°R) 

=  8.315  kJ/(kg-moI)  (K) 

-  1545  ft-lbf  (lb-mol)  (8R) 

=  1.986  Btu/(lb-mol)  UR) 
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SECTION  I 


INTRODUCTION 

Two  of  the  primary  requirements  of  a  gas  turbine  combustor  design  are  to  maintain 
combustion  over  a  wide  range  of  operating  conditions  (including  the  transient  state  of  rapid 
deceleration)  and  to  achieve  rapid  relighting  and  engine  acceleration  after  flame-out  at  altitude. 
The  modern  military  fighter  and  commercial  engines  are  called  upon  to  operate  at  very  low  inlet 
temperatures,  pressures,  and  at  primary  zone  fuel  air  ratios  that  on  the  average  are  well  outside 
the  limits  of  flamability  for  hydrocarbon  air  mixtures.  Under  these  conditions,  combustion  may 
become  unstable  which  results  in  a  tendency  for  the  flame  to  blow-out  and  not  relight. 

Present  and  future  generations  of  gas  turbine  combustors  are  being  designed  for  higher 
turbine  inlet  temperature  and  minimum  smoke  emission.  This  requires  that  a  larger  percentage 
of  the  air  be  routed  to  the  combustor  primary  zone  to  obtain  proper  stoichiometry  during  high 
power  operation.  However,  during  altitude  ignition  and  engine  deceleration,  the  primary  zone 
stoichiometry  becomes  excessively  lean.  This  reduces  combustion  stability  margin  and  leads  to 
potential  deceleration  blowout  and  an  inability  to  ignite  at  required  altitude. 

A  number  of  experimental  evaluations  to  determine  altitude  ignition  and  lean  deceleration 
blowout  limits  have  been  sponsored  by  the  Air  Force,  Navy,  NASA  and  Pratt  &  Whitney  over  the 
past  few  years.  These  tests  have  provided  a  wealth  of  data.  The  data  sources  selected  for  this 
study  are  as  shown  in  Table  2.  The  objective  of  the  Altitude  Ignition  and  Lean  Decel  Study  has 
been  to  use  this  data  base  and  to  develop  analytical  tools  consisting  of  new  or  improved  models 
for  predicting  altitude  ignition  and  lean  decel  performance.  The  models  may  then  be  used  as  the 
basis  for:  (1)  serving  as  a  design  tool  for  new  combustion  systems;  (2)  permitting  necessary 
design  modifications  to  be  made  in  the  early  stages  of  the  combustor  design;  and  (3)  reducing  the 
time  and  cost  for  testing  new  combustion  systems  in  order  to  define  ignition  and  stabilization 
limits. 

The  Altitude  Ignition  and  Lean  Decel  Study  was  conducted  in  three  tasks  as  listed  below. 

•  Task  I  —  An  in-depth  literature  survey  was  made  of  existing  ignition  and 
stability  theory  emphasizing  altitude  ignition  and  lean  decelera¬ 
tion  blowout 


•  Task  II  —  The  available  in-house  and  external  engine  test  and  rig  test 

data  were  gathered  to  support  the  analytical  effort.  Emphasis 
was  placed  on  engine  data  because  it  includes  effects  such  as 
combustor  inlet  air  flow  distortion  and  control  acceleration  and 
deceleration  schedules,  both  of  which  are  difficult  to  duplicate 
in  a  rig 

•  Task  III  —  The  selection  and  development  of  the  model  which  best  fits 

the  data  for  altitude  ignition  (windmill  and  spooldown  modes) 
and  the  best  correlation  for  lean  deceleration  blowout  was 
made. 


mmimsm 
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TABLE  2. 


ALTITUDE  IGNITION  AND  LEAN  DECEL  DATA 


Engine 

Type  of  Data  Gathered 

a.  Military  Engines 

F100(3)  Engine  P1068X-13 

Altitude  Ignition  (Windmill  and  Spooldown) 

F  100(3)  Engine  FX2 19-21 

Altitude  Ignition  (Spooldown) 

PW1128  Development  Engine  FX225-19 

Altitude  Ignition  (Windmill  and  Spooldown) 

ILC  Core  Engine  FX208-29 

Altitude  Ignition  (Windmill) 

TF30- P-4 14/4 14A  Engine  P-286-1  and  2 

Altitude  Ignition  (Windmill) 

FI00-PW-200  Engine  FX217-19 

Lean  Decel  Blowout 

F100  ILC  Development  Engine  FX217-21 

lean  Decel  Blowout 

b.  Commercial  Engines 

PW2037  Engine  X-669-1/B747 

Altitude  Ignition  (Windmill) 

JT9D-7R4  Engine  X-579-31 

Altitude  Ignition  (Windmill) 

JT9D-7R4  Engine  X-491-2 

Altitude  Ignition  (Windmill) 

PW2037  Engine  X-663-1 

Lean  Decel  Blowout 

PW2037  Engine  X-664-7 

Lean  Decel  Blowout 

I4U6T 


The  development  of  altitude  ignition  and  lean  deceleration  under  Task  III  was  conducted 
according  to  the  following  steps.  First,  relationships  were  established  which  associate  character¬ 
istic  times,  with  the  physical  processes  in  the  combustor  for  ignition  and  stabilization.  Second, 
these  expressions  were  interested  according  to  the  criteria  for  successful  ignition  and 
stabilization.  Third,  statistical  analysis  was  used  to  compare  the  various  expressions  and  select 
the  appropriate  terms  that  formulate  models  which  best  fit  the  data  for  altitude  ignition  and  lean 
deceleration  blowout.  These  models  have  been  based  entirely  on  full-scale  engine  tests. 
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SECTION  II 

DEVELOPMENT  OF  SP/tRK  IGNITION  AND  FLAME  STABILIZATION  MODELS 


1.  DEFINITION  OF  PARAMETERS 

The  success  in  adapting  existing  ignition  and  stability  models  to  engine  operating 
conditions  was  dependent  in  a  large  part  to  properly  defining  the  various  fuel  and  air  properties, 
operating  and  design  parameters  which  were  employed.  The  definition  of  these  parameters  is 
presented  in  this  section.  In  most  cases  the  basis  for  defining  the  parameter  was  clear  cut.  In 
other  cases  the  definition  was  obtained  through  a  systematic  evaluation  of  the  effect  of  the 
parameters  involved. 

a.  Fuel  Properties 

For  the  military  and  commercial  engine  data  in  this  study,  expressions  were  developed  for 
the  following  fuel  properties: 

•  Density 

•  Surface  tension 

•  Absolute  viscosity 

•  Vapor  pressure 

•  Latent  heat  of  vaporization 

•  Specific  heat 

•  Mass  transfer  number. 

The  fuel  properties  were  determined  as  a  function  of  the  fuel  temperature.  The  fuel 
temperature  was  obtained  from  the  data  set  plus  an  estimated  temperature  rise  in  the  fuel  system 
and  injector.  In  addition,  some  fuel  properties  were  evaluated  as  a  function  of  10  percent 
distillation  temperature  for  spark  ignition  phenomena  and  50  percent  distillation  temperature 
for  flame  stabilization  phenomena. 

b.  Air  and  Gas  Properties 

The  air  and  gas  properties,  such  as  density,  absolute  viscosity,  specific  heat,  thermal 
diffusivity,  thermal  conductivity  and  Prandtl  number  were  determined  as  a  function  of 
combustor  inlet  temperature  or  gas  temperature  in  the  primary  zone. 

c.  Injector  Velocity  and  Spray  Angle 

The  injector  velocity  was  evaluated  by  considering  the  fuel  injector  and  swirler  flow  as  one 
dimensional  incompressible  flow  with  axial  and  tangential  components.  Given  the  injector  and 
swirler  air  flow  rate,  the  effective  and  physical  flow  areas,  combustor  pressure  and  temperature, 
it  is  possible  to  determine  the  average  axial  component  of  the  injected  velocity  and  the  spray 
angle  of  the  injector.  The  axial  velocity  and  the  spray  angle  were  used  to  calculate  the  tangential 
component,  and  the  addition  of  these  components  provided  on  average  injected  velocity  at  the 
exit  plane  of  the  injector.  Also,  reference  velocities  were  calculated  at  the  plane  of  the  ignitor  and 
primary  holes. 

d.  Sauter  Mean  Diameter 

To  calculate  the  droplet  evaporation  time  for  a  particular  combustor,  the  fuel  spray  droplet 
Sauter  mean  diameter  (SMD)  must  be  determined  for  the  particular  combustion  system  of 
interest.  Unfortunately,  droplet  size  measurements  were  not  available.  Instead,  empirical 
correlation  parameter  were  used  to  predict  the  SMD  for  each  atomizing  nozzle. 


3 


(4*71 


(1)  Airblast  Atomizing  Nozzles 


The  F100(3),  PW1128,  and  PW2037  engines  employ  airblast  atomizing  fuel  nozzles. 
Rizkalla  and  Lefebvre’s  correlation  parameter  for  this  type  of  atomizer,  Reference  30,  was  used 
to  estimate  the  relative  fuel  spray  droplet  diameter. 

«■> -  Kfcf  (*■)“(>♦  ^;)>]  ♦  [n,©'V^)c]  (1) 

where: 

N,,  N2  =  Constant  (dependent  on  fuel  nozzle  design) 
pa  =  Air  density,  (kg/m3) 

pf  =  Fuel  density,  (kg/m3) 

of  =  Surface  tension,  (N/m) 

pf  =  Dynamic  viscosity,  (kg/m-s) 

Wf  =  Fuel  flow,  (kg/s) 

Wa  =  Air  flow,  (kg/s) 

U;nj  =  Injected  velocity,  (m/s) 

tp  =  Liquid  film  thickness  (microns). 

(2)  Hybrid  Atomizing  Nozzles 

The  JT9D-7R4  combustor  employs  hybrid  atomizing  fuel  nozzles.  For  this  type  of  atomizer, 
a  correlation  parameter  developed  at  Pratt  &  Whitney  was  used  to  predict  the  relative  fuel  spray 
droplet  diameter. 

SMD  =  [  N,  (W rf,R  [e*  "“^'f^j/fAPff442 


where: 


N3  =  Constant  (dependent  on  fuel  nozzle  design) 

Wf  =  Fuel  flow  (kg/s) 

Tf  =  Fuel  temperature  (°K) 

APf  =  Pressure  drop  in  the  fuel  nozzle  passage  (kPa). 

(3)  Pressure  Atomizing  Nozzles 

The  TF30  combustor  employs  pressure  atomizing  fuel  nozzles.  Lefebvre  correlation 
parameter  for  this  type  of  atomizer  (Reference  31)  was  used  to  estimate  the  relative  fuel  spray 
droplet  diameter. 


SMD  -  N2  o?25  vf*  Wfv2r’  APfn'5  p„0  2R 


where: 


N2  =  Constant  (dependent  on  fuel  nozzle  design) 
vf  «  Kinematic  viscosity,  (m2/sec) 
p  =  Density  of  the  gas,  (kg/m3). 


4*7C 


e.  Effective  Fuel  Air  Ratio  and  Gas  Temperature  in  Primary  Zone 


The  rate  of  fuel  evaporation  is  not  always  sufficiently  high  to  ensure  that  all  of  the  fuel  is 
fully  vaporized  within  the  primary  zone.  Thus,  an  iterative  procedure  was  developed  to  determine 
the  fraction  of  fuel  vaporized  or  an  “effective”  fuel  air  ratio  and  the  corresponding  gas 
temperature  within  the  combustion  zone.  The  procedure  starts  with  an  initial  guess  effective  fuel 
air  ratio  equal  to  the  stoichiometric  fuel  air  ratio  and  determines  the  corresponding  gas 
temperature.  The  fraction  of  fuel  vaporized  is  then  calculated  according  to  the  following 
expression  (Reference  31): 


f-  -  8  (£)  (i).  ■"(>  -  B)  (^T^r)  (1  -  0.22R4’) 


where: 


(4) 


if 


ft1 

pz 

SMD 

Ren 


Fraction  of  fuel  vaporized  within  combustion  zone 
Gaseous  thermal  conductivity,  (J/m-s-K) 

Gaseous  specific  heat,  (J/kg-K) 

Mass  transfer  number  based  on  50  percent  boiling  temperature 
Volume  of  the  primary  zone,  m3 
Primary  zone  airflow,  (kg/sec) 

Sauter  mean  diameter,  microns 

Droplet  Reynolds  number  based  on  mean  drop  size 


From  the  fraction  of  fuel  vaporized,  an  effective  fuel  air  ratio  is  calculated  and  compared  to 
the  initial  effective  fuel  air  ratio.  If  the  calculated  value  of  effective  fuel  air  ratio  is  not  equal  to 
the  initial  value  used  to  obtain  the  gas  temperature,  an  iteration  process  is  used  to  balance  them. 
It  should  be  noted  that  equation  4  allows  the  fraction  of  fuel  vaporized  to  exceed  unity  which 
means  that  the  fuel  is  fully  vaporized  within  the  recirculation  zone.  In  these  circumstances,  the 
primary  zone  fuel  air  ratio  is  used  to  determine  the  gas  temperature  and  properties  in  the 
primary  zone. 


f.  Primary  Zone  Volume 

The  primary  zone  volume  is  required  to  evaluate  the  fraction  of  fuel  vaporized  and  the 
volumetric  hot  residence  time.  The  volume  of  the  primary  zone  of  an  annular  combustor  formed 
by  the  spray  angle  of  the  fuel  injector  and  the  primary  zone  jets  is  calculated  in  the  following 

manner.  The  volume  is  assumed  to  take  the  shape  of  a  frustum  of  a  right  circular  cone  with  the 

slant  height  at  the  spray  angle,  the  lower  base  at  the  plane  of  the  primary  holes  and  the  upper 
base  at  the  exit  of  the  injector  plane.  It  uses  the  primary  zone  length  as  the  vertical  height  of  the 
frustum.  While  the  radius  of  the  lower  and  upper  bases  are  represented  by  one  half  of  the  annulus 
height  in  the  plane  of  the  primary  holes  and  the  exit  swirler  radius  of  the  injector,  respectively. 

g.  Spark  Kernel  Quench  Diameter 

A  definition  of  minimum  ignition  energy  is  required  to  evaluate  the  diameter  of  the  spark 
kernel  (dq).  As  suggested  by  Ballal  and  Lefebvre,  Reference  1,  minimum  ignition  energy  is 
defined  as  the  energy  required  to  heat  a  spherical  volume  of  air  with  diameter  dq,  up  to  the 
stoichiometric  flame  temperature,  which  results  in  the  following  equation: 

.  r  Em,n  r 

q  L  7i  6  p,  c,.„  AT„  J  (fl) 

where: 

Kmjn  =  Minimum  ignition  energy,  (J) 
cpi)  “  Specific  heat  of  air,  (J/kg-K) 

ATS,  =  Stoichiometric  (lame  temperature  rise,  (K) 
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h.  Evaporation  Coefficient 


To  compute  the  droplet  evaporation  time  in  the  spark  ignition  and  flame  stabilization 
models  developed  by  Peters,  Plee,  and  Mellor,  Reference  2,  it  is  necessary  to  determine  the 
droplet  evaporation  coefficient  which  accounts  for  fuel  properties  influences  and  is  corrected  for 
convection  according  to  the  phenomena  of  interest. 

For  the  ignition  model 
8k 

Pi  =  In  (1+B)  (0.185  Re°D6) 

fiS.d  (6) 


where: 


p.g 

Pf 

B 

Re 


D 


Gaseous  thermal  conductivity  (J/m-s-K) 

Gaseous  specific  heat  (J/kg-K) 

Fuel  density  at  10  percent  boiling  point  kg/m3 

Transfer  number  based  on  the  10  percent  boiling  temperature 

Droplet  Reynolds  number  based  on  mean  drop  size. 


For  the  stabilization  model 


Ps  = 


8k„ 


PfCp 


In  (1+B)  (1+0.276  ReJ,'*  Sc,J) 


(7) 


where: 


Pf 

B 

Sc 


=  Fuel  density  at  50  percent  boiling  point,  (kg/m3) 

-  Transfer  number  based  on  50  percent  boiling  point  temperature 

-  Gas  phase  Schmidt  number. 


Root  Mean  Square  of  the  Fluctuating  Velocity 


Turbulence  measurements  were  not  available  from  the  data  presented  in  Section  I. 
However,  using  the  approach  of  Lefebvre,  Reference  28,  the  root  mean  square  of  the  fluctuating 
velocity,  u',  was  predicted  from  the  pressure  loss  across  the  particular  combustor  and  the  inlet 
temperature  and  pressure.  The  pressure  energy  loss  must  appear  as  turbulence  energy  prior  to 
dissipation  as  molecular  motion.  If  the  pressure  loss  is  proportional  to  the  turbulent  kinetic 
energy.  In  equation  form: 


AP  cc  -g-  p,  (u'f 


(8) 


Substituting  the  perfect  gas  law  equation  for  the  density  and  solving  for  u'  equation  (8) 
becomes: 


u'  =  e[l.5RuT1((AP/P)/100)J 


(9) 


where: 


e 

Ru 

Ma 

T, 

AP 

V 


Constant  of  proportionality 
Universal  gas  constant  8815  (.I/kg  mole  K) 
Molar  mass  of  air,  29  (kg/kgmole) 
Combustor  inlet  temperature  (K) 


Percent  combustor  pressure  drop. 
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j.  Laminar  and  Turbulent  Flame  Speeds 


The  laminar  flame  speed  is  calculated  from  a  correlation  developed  by  Pratt  &  Whitney. 
This  correlation  was  calibrated  against  ethane  and  propane  data,  and  includes  the  effects  of 
pressure,  inlet  temperature  and  fuel  air  ratio. 


S,. 


exp  (-2xl04/Tad  +2xl04/TadJ 


where: 


SL  =  Laminar  flame  speed,  (m/s) 

T.,  =  Combustor  inlet  temperature  (R) 

Tad  =  Adiabatic  flame  temperature  (R) 

o  =  Reference  conditions. 


(10) 


The  turbulent  flame  speed  is  based  on  the  Shchelkim's  model  which  relates  turbulence 
properties  to  the  increase  in  specific  surface  area  of  the  flame.  This  approach  leads  to  a 
relationship  of  the  form: 


where: 


ST=  Turbulent  flame  speed,  (m/s) 
in  which  q  is  a  constant  of  the  order  of  unity. 

2.  LITERATURE  SURVEY 

An  indepth  literature  survey  was  conducted  of  existing  ignition  and  stability  theory  and 
correlations  emphasizing  altitude  ignition  and  lean  deceleration  blowout.  Titles,  authors  sources, 
and  a  brief  description  of  literature  reports  applicable  to  this  study  are  shown  in  Appendix  A.  A 
list  of  all  available  reports  uncovered  during  the  literature  survey  can  be  found  in  the  Reference 
Section.  The  following  models  were  considered  in  this  study. 

3.  SPARK  IGNITION  MODELS 

The  spark  ignition  model  developed  by  Peters  and  Mellor,  Reference  2,  states  that  for 
ignition  to  occur,  the  energy  of  the  spark  must  heat  an  initial  volume  such  that  the  heat  release 
rate  within  that  volume  is  greater  than  the  loss  rate.  Heat  generation  is  limited  first  by  a  droplet 
evaporation  rate  and  then  by  a  chemical  reaction  rate.  The  heat  loss  in  a  gas  turbine  combustor 
environment  is  due  to  turbulent  diffusion;  hence,  it  is  controlled  by  a  turbulent  mixing  time.  The 
ignition  limit  is  reached  when  ts1,  the  turbulent  mixing  time,  equals  the  sum  of  teb,  the  droplet 
evaporation  time  plus  the,  the  chemical  reaction  time.  In  equation  form: 


t«|  the  teh 


(12) 


The  spark  ignition  model  by  Ballal  and  Lefebvre,  Reference  1,  states  that  for  ignition  to 
occur,  the  passage  of  the  spark  must  create  a  small,  roughly  spherical  volume  of  air  (spark  kernel) 
whose  temperature  is  sufficiently  high  to  initiate  evaporation  and  chemical  reaction.  If  the  rate 


of  heat  release  by  combustion  exceeds  the  rate  of  heat  loss  by  thermal  conduction  and  turbulent 
diffusion,  then  the  kerne!  grows  to  a  size  capable  of  flame  propagation.  The  ignition  limit  is 
reached  when  the  time  required  for  the  fuel  to  evaporate  (xei)  chemically  react  (xcl)  is  equal  to  or 
less  than  the  time  required  for  the  cold  mixture  to  quench  (xq)  the  spark  kernel  by  thermal 
conduction  and  turbulent  diffusion.  In  equation  form: 


T  a  T  +T 

**q  *,ei  1  *’Cl 


a.  Evaluation  of  the  Characteristic  Times 


(1)  Turbulent  Mixing,  xs/ 


The  turbulent  mixing  time,  xg!,  was  defined  as  the  ratio  of  the  diameter  of  the  spark  kernel 
(dq)  to  the  mean  flow  velocity  at  the  plane  of  the  ignitor  (Uign).  In  equation  form: 

T-si  -  dq/'Ujp,  (msec)  (14) 

(2)  Spark  Kernel  Quench  Time,  xq 

The  spark  kernel  quench  time,  xq ,  was  defined  as  the  ratio  of  the  heat  capacity  of  the  spark 
kernel  divided  by  the  average  rate  of  heat  loss  from  the  kernel  by  thermal  conduction  and 
turbulent  diffusion,  i.e., 

For  moderate  turbulence 


).68u'3LT(msec) 


For  high  turbulence,  this  reduces  to: 


T«  “  0.6W'(msec> 


where: 

dq  =  Spark  kernel  diameter,  (m) 

a  =  Thermal  diffusivity,  (m2/s) 

u'  =  Root  mean  square  of  the  fluctuating  velocity,  (m/s) 

(3)  Fuel  Evaporation  Times,  xrt  and  x^ 

Peters  and  Mellor  computed  the  droplet  evaporation  time,  xeb,  from  the  “d2  law"  based 
upon  the  initial  SMD  of  the  spray.  The  droplet  evaporation  time  was  defined  as  the  ratio  of  the 
square  of  the  Sauter  mean  diameter  divided  by  the  droplet  evaporation  coefficient.  Since  the 
total  amount  of  fuel  being  vaporized  contributes  to  the  ignition  process,  xeh  was  divided  by  the 
total  number  of  drops  in  the  spark  kernel.  Peters  and  Mellor  noted  that  the  equivalence  ratio  in 
the  primary  zone  is  proportional  to  the  total  drops  in  the  spray.  In  equation  form: 


x,b  =  (SMD)7Pj  q>ien  (msec) 


14*7< 


where: 


SMD  =  Sauter  mean  diameter  of  the  spray,  (m) 

Pi  =  Droplet  evaporation  coefficient  for  spark  ignition,  (m2/s) 

<pign  =  Equivalence  ratio  near  the  ignitor. 

Ballal  and  Lefebvre  defined  the  evaporation  time,  iei,  as  the  difference  between  the  fuel 
liquid  evaporation  time  and  the  time  to  produce  the  vapor  that  already  exists  in  the  system.  Fuel 
liquid  evaporation  time  was  obtained  as  the  ratio  of  the  mass  of  fuel  contained  within  the  spark 
kernel  to  the  rate  of  fuel  evaporation,  while  vaporization  time  was  obtained  as  the  ratio  of  the 
mass  of  vapor  present  to  the  rate  of  evaporation  of  the  liquid  drops  in  the  spray.  Therefore, 
evaporation  time  was  calculated  in  the  following  manner: 


For  high  turbulence  conditions 


(Pr/2)(1  -  ft)  pf  (SMD)1 5 
Z  P»S  Ma  S  U7"  ’  (Pm,  In  (1  +  B„] 


■  [msec) 


where: 


Prandtl  number 

Fraction  of  total  fuel  in  form  of  vapor 
Fuel  density,  (kg/m3) 

Sauter  mean  diameter,  (m) 

C,  c2°"Vc3:i 

Ratio  of  the  mean  surface  area  diameter  to  the  Sauter  mean  diameter 
Ratio  of  the  mean  diameter  to  the  Sauter  mean  diameter 
Ratio  of  the  mean  volume  diameter  to  the  Sauter  mean  diameter 
Air  density  (kg/m:i) 

Dynamic  viscosity,  (kg/m-sec) 

Root  mean  square  of  the  fluctuating  velocity,  (m/sec) 

Mass  transfer  number  at  stoichiometric  temperature 
Injector  equivalence  ratio. 


(4)  Chemical  Reaction  Times,  xhc  and  xcl 


Peters  and  Mellor  derived  the  chemical  reaction  time,  xhf  from  the  Arhenius  expression  in 
terms  of  the  rate  of  hydrocarbon  oxidation. 


A  exp  (E  RT„ 


-  (msec) 


where: 

A  =  Pre-exponential  factor,  10  s  (msec-kg/m3) 

E  =  Activation  energy,  26,100  (cal/mole) 

R  =  Universal  gas  constant,  1.986  (cal/mole- K) 

Tst  =  Stoichiometric  adiabatic  flame  temperature,  (K) 
pK  =  Gas  phase  density,  (kg/m3) 

‘Pi^n  =  Equivalence  ratio  near  the  igniter. 

Ballal  and  Lefebvre  described  the  chemical  reaction  time  as  a  function  of  the  laminar  and 
turbulent  flame  speeds. 


For  high  turbulence  conditions: 

15.6a  ,  . 

T"  "  u'(ST  -  0.63  u'T  (msec>  (19) 

where: 

a  =  Thermal  diffusivity,  (m2/s) 

u'  =  Root  mean  square  of  the  fluctuating  velocity,  (m/s) 

S-p  =  Turbulent  flame  speed,  (m/s). 

4.  FLAME  STABILIZATION  MODELS 

In  the  present  study,  the  flame  stabilization  process  is  viewed  as  occurring  in  the  shear 
layer  surrounding  the  combustor  recirculation  zone  by  way  of  turbulent  mixing  of  fresh  reactants 
and  hot  products  and  partially  oxidized  fuel.  The  stability  limit  is  reached  when  the  volumetric 
residence  time,  xvhr,  of  the  hot  turbulent  eddies  present  in  the  shear-layer  region  is  of  the  order  of 
evaporation  time,  xes,  plus  chemical  reaction  time,  xC8.  In  equation  form: 

*vhr  =  tfS  +  x„  (msec) 

(20) 

Plee  and  Mellor  considered  the  same  approach  and  criteria  with  the  exception  that  the 
volumetric  hot  residence  time  was  replaced  with  the  turbulent  mixing  time.  Also,  the  turbulent 
mixing  time  was  multiplied  by  the  ratio  of  T3/Tst.  This  accounts  for  the  fluid  mechanic 
acceleration  across  the  shear  layer  due  to  the  temperature  gradient,  and  was  included  in  the 
characteristic  times  of  evaporation,  xEB,  and  chemical  reaction  xHC  so  that  the  fluid  mechanics 
and  chemistry  are  uncoupled.  In  equation  form: 

tsi.  =  Tkb  +  *HC  (msec) 

(21) 

a.  Evaluation  of  the  Characteristic  Times 


(1)  Volumetric  Hot  Residence  Time,  xyhr 

The  volumetric  hot  residence  time,  xvhr,  is  defined  as  the  ratio  of  the  primary  zone  volume 
divided  by  the  airflow  through  the  primary  zone  times  the  gas  density.  In  equation  form: 

*vhr  =  Vp2pK/Wap2(msec) 


where: 

Vpz  =  Volume  of  the  primary  zone,  (m3) 
p  =  Gas  density,  (kg/m3) 

Wapz  =  Air  in  the  primary  zone,  (kg/sec). 
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(2)  Turbulent  Mixing  Time,  xSL 

The  turbulent  mixing  time,  tsl  was  characterized  by  an  eddy  decay  time  which  was  related 
to  the  macroscopic  characteristic  length  for  CO  emissions,  lr0,  and  reference  velocity  in  the 
primary  zone,  U  ,  introduced  by  Mellor  and  Washam  to  model  emissions: 

tsi.  =  tsl =  I,-,,/ U,«  (msec) 


The  characteristic  length  for  CO  emissions  is  defined  by: 


lc<>  =  [d,,!„,b  +  1,,,'j  ' 


where: 


1  -•  Axial  distance  from  the  discharge  plane  of  the  fuel  injector  to 

the  center  line  of  the  primary  holes 

d(.l(mt)  =  Combustor  diameter. 

(3)  Fuel  Evaporation  Times,  tes  and  xES 

The  droplet  evaporation  time,  tes  is  defined  as  the  ratio  of  the  mass  of  fuel  present  to  the 
average  rate  of  fuel  evaporation  for  a  polydisperse  spray.  Thus  according  to  Ballal  and  Lefebvre, 
the  droplet  evaporation  time  is  given: 


Ci!  Pr  cp,B  (SMDf  (l-£2) 


8C,  kK  <p,.  In  (1+B)  (1+0.25  Cf  Re'/,')  ' 


(msec) 


where: 


pf  =  Fuel  density  at  50  percent  distillation  temperature,  (kg/m!) 

c  =  Specific  heat  of  gas,  (J/kg-K) 

SMD  =  Sauter  mean  diameter,  (m) 

S2  =  Fraction  of  total  fuel  in  form  of  vapor 

C,  =  Ratio  of  the  mean  surface  area  diameter  to  the  Sauter  mean 

diameter 

C2  =  Ratio  of  the  mean  diameter  to  the  Sauter  mean  diameter 

C;!  =  Ratio  of  the  mean  volume  diameter  to  the  Sauter  mean  diameter 

kK  =  Thermal  conductivity  of  gas,  (J/m-s-K) 

*Ppz  ~  Equivalence  ratio  in  the  primary  zone 

B  =  Mass  transfer  number  at  50  percent  distillation  temperature 
ReD  -  Droplet  Reynolds  number  based  on  mean  drop  size. 

Plee  and  Mellor  computed  the  droplet  evaporation  time,  using  equation  (16).  However, 
the  convection  correction  in  the  evaporation  coefficient  was  modified  according  to  equation  (7), 
the  equivalence  ratio  was  removed  and  properties  were  evaluated  at  50  percent  distillation 

temperature  for  the  phenomena  of  flame  stabilization.  This  results  in  the  following  equation 

form: 


(SMDf  ().  (  qf'-)  l msec) 


■'  »  ■*»  •* « "?•  i .“T -  \  -*■  ■a.  *■“«  •*"„  •»".  \  •* -  - v  «■„  v,  *  «*,  •-  «~N  - 


where: 


SMD  =  Sauter  mean  diameter,  (microns) 

Ps  =  Evaporation  coefficient  for  stability,  (m2/s) 

Tst  =  Stoichiometric  adiabatic  flame  temperature,  (K) 

T3  =  Combustor  inlet  temperature,  (K). 

(4)  Chemical  Reaction  Time,  tc#  and  thc 

The  chemical  reaction  time,  tC9,  considered  in  this  study  for  the  phenomena  of  flame 
stabilization  is  defined  as  in  the  spark  ignition  phenomena  developed  by  Ballal  and  Lefebvre 
equation  (19).  Flame  stabilization  and  ignition  are  similar  phenomena.  In  both  cases,  a  fuel  and 
air  mixture  must  be  heated  so  that  the  fuel  evaporates,  mixes  with  the  air  and  chemical  reactions 
begin  at  the  rate  sufficient  for  establishing  flame.  The  difference  between  ignition  and  flame 
stabilization  is  the  energy  source  that  utilizes  combustion,  hot  recirculation  gas  for  flame 
stabilization,  and  a  spark  ignition.  Consequently,  the  formulation  of  the  chemical  reaction  time 
for  ignition  is  the  same  for  the  flame  stabilization. 

The  chemical  reaction  time,  thc,  developed  by  Plee  and  Mellor  for  the  phenomena  of  flame 
stabilization  is  defined  in  the  same  manner  as  in  the  spark  ignition  phenomena  with  the 
exception  of  the  pre-exponential  factor  is  10-4 ,  the  activation  energy  is  21,000  cal/mole,  the  gas 
density  term  is  removed  and  the  equivalence  ratio  is  evaluated  as  the  primary  zone  equivalence 
ratio.  In  equation  form: 

tH(-  =  (A/«ppi)  exp  (E/RUTJ  (-^-)  (msec) 

where: 

A  =  Pre-exponential  factor,  10-4  (msec-kg/m3) 

E  =  Activation  energy,  21,000,  (cal/mole) 

<ppz  =  Equivalence  ratio  in  the  primary  zone 
Ru  =  Universal  gas  constant,  1.986,  (cal/mole-K) 

Tst  =  Stoichiometric  adiabatic  flame  temperature,  (K) 

T3  =  Combustor  inlet  temperature,  (K). 


SECTION  III 


ALTITUDE  IGNITION  MODELS 


Altitude  ignition  followed  by  engine  acceleration  has  been  one  of  the  most  significant 
combustor  performance  variables  affected  by  the  trend  toward  increased  air  flow  in  the  primary 
zone.  In  recent  years,  analytical  and  experimental  studies  have  been  conducted  to  provide  a 
useful  theoretical  foundation  for  relating  ignition  characteristics  to  combustor  operating 
parameters.  Some  of  these  theories  were  presented  in  Section  II,  and  were  considered  in  the 
development  of  altitude  ignition  models  for  windmill  and  spooldown  engine  performance. 

Altitude  ignition  models  were  developed  using  F100(3),  PW1128  Development,  F100-ILC 
Core,  PW2037,  JT9D-7R4  and  TF30-P-414/414A  engine  data.  Ignition  operating  conditions 
were  determined  by  plotting  the  high  pressure  rotor  speed  and  combustor  pressure  as  a  function 
of  time.  Ignition  followed  by  engine  acceleration  was  identified  by  the  sudden  increase  in 
combustor  pressure  and  verified  by  the  increase  in  high  pressure  rotor  speed.  However,  if  no 
response  occurred,  the  operating  point  was  defined  as  no  ignition.  Also,  blowout  points  were 
determined  by  the  sudden  decrease  in  combustor  pressure  or  high  pressure  rotor  speed  during 
engine  acceleration. 

The  process  of  altitude  ignition  is  viewed  to  occur  in  the  following  manner.  Upon 
attempting  to  relight  an  engine  either  by  windmill  or  spooldown  mode,  the  spark  kernel  must 
grow  to  a  minimum  spherical  volume  of  sufficiently  high  temperature  to  initiate  rapid 
evaporation,  and  sustain  both  chemical  reaction  and  flame  propagation  while  the  primary  zone 
operating  conditions  must  be  stable.  Otherwise,  ignition  followed  by  engine  acceleration  will  not 
occur. 


Comparison  of  the  degree  of  correlation  of  the  ignition  and  flame  stabilization  models  is 
presented  in  Section  II.  The  characteristic  times  developed  by  Ballal  and  Lefebvre  for  the 
phenomena  of  spark  ignition  were  selected  to  describe  the  physical  processes  in  the  combustor 
during  windmill  and  spooldown  ignition.  The  spark  kernel  ignition  model  includes  the  effects  of 
thermal  and  turbulent  diffusion,  evaporation  and  chemical  reaction  time.  The  following 
expressions  describe  the  characteristic  times  associated  with  each  physical  process: 


Spark  kernel  quenching  time: 
dq 

u  -  a64Tr(msec) 

Evaporation  time: 

_  Cti(Pr,2)(l-fl)pf  (SMD)15 
T"  ~  Z  p“  p"  s  u'"s  <p,n)  In  (1+BJ 


(msec) 


Chemical  reaction  time: 


Cci  (15.6a)  , 

Tli  “  u'  (ST-0.63  u'I  (msec) 


(27) 


(28) 


(29) 


Where  the  proportionality  and  the  constant  weighting  factors,  Cei  =  Crj  =  2.66  and  1.0  for 
windmill  and  spooldown  ignition,  are  necessary  because  characteristic  times  are  simply  order  of 
magnitude  estimates  of  the  processes  involved. 
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Since  successful  ignition  must  be  followed  by  engine  acceleration,  the  primary  zone  stability 
limit  was  also  analyzed  at  the  time  of  ignition.  The  flame  stabilization  and  spark  ignition  are 
similar  phenomena.  Consequently  the  formulation  of  the  physical  processes  for  the  flame 
stabilization  model  is  very  similar  to  the  spark  ignition  model.  The  primary  zone  flame 
stabilization  model  considers  the  effects  of  volumetric  flow  residence,  evaporation  and  chemical 
reaction.  These  physical  processes  are  defined  analytically  as: 


Volumetric  hot  residence  time: 


T 


vhr 


v,»  P, 


(msec) 


(30) 


Evaporation  time: 

T  =  C,C;PfCM(SMDf(l-fl) 
ps  8  C,  kK<ppz  In  (1+0.25  C“-25  Re?,*)  (mSeC) 


Chemical  reaction  time: 


C„(15.6a) 

T"  "  u7^  -  0.63  u'f  (msec) 


where,  C^  =  Ccs  =  1.66  for  windmill  and  spooldown. 

1.  WINDMILL  IGNITION  AND  PRIMARY  ZONE  FLAME  STABILIZATION 


(31) 


(32) 


In  the  present  study,  windmill  ignition  limit  data  from  seven  engines  F100(3),  PW1128 
Development,  F100-ILC  Core,  PW2037,  JT9D-7R4  and  TF30-P-414/414A  were  analyzed.  Linear 
correlation  models  of  windmill  ignition  and  primary  zone  stability  limits  were  obtained  as  shown 
in  Figures  1  and  2. 


A  good  correlation  model  of  the  spark  ignition  limit  was  obtained  by  correlating  the 
available  spark  kernel  quench  time  tq  as  a  function  of  the  evaporation  plus  chemical  reaction 
time  tei  +  xci.  The  solid  line  in  Figure  1  represents  the  spark  ignition  limit  and  separates  the 
regions  of  ignition  and  no  ignition.  With  increased  evaporation  plus  chemical  reaction  time  for 
specified  spark  quench  time,  the  engine  approaches  the  limit  of  no  ignition,  and  as  a  result, 
longer  spark  quench  time  is  required  to  establish  successful  ignition. 

Figure  2  presents  the  relationship  between  volumetric  hot  residence  time  and  the  sum  of 
evaporation  plus  chemical  reaction  for  the  phenomena  of  primary  zone  flame  stabilization 
immediately  after  ignition.  The  solid  line  represents  the  primary  zone  stability  limit  and 
separates  the  regions  of  stable  and  unstable  conditions.  Current  and  advanced  engines  with 
increased  primary  zone  air  flow  and/or  reduced  volume  approach  the  limit  of  stability  by 
decreasing  volumetric  hot  residence  time  for  a  specified  chemical  plus  evaporation  time.  To 
obtain  successful  ignition  followed  by  engine  acceleration,  the  combustor  operating  conditions 
must  exist  within  the  ignition  and  stable  region. 
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The  equations  of  the  limit  lines,  coefficients  of  determination  (0.920  and  0.950)  and 
standard  deviations  (3.725  and  0.037)  for  ignition  and  stability  are  shown  in  Figures  3  and  4, 
respectively.  The  dotted  lines  represent  the  95  percent  confidence  band.  For  any  value  of  the 
evaporation  plus  chemical  reaction,  the  interval  about  the  Y  axis  between  the  dotted  lines 
represents  the  true  value  of  Y  with  95  percent  confidence  band. 

The  correlation  models  are  valid  over  a  wide  range  of  combustor  operating  conditions  which 
include  pressures  from  0.2  to  1.0  atm,  primary  zone  velocity  of  5  to  40  ft/sec,  Sauter  mean 
diameter  of  30  through  180  microns,  primary  zone  equivalence  ratio  of  0.4  to  2.8. 

The  evaporation  and  chemical  reaction  effects  on  windmill  spark  ignition  limit  are  shown 
in  Figures  5  and  6,  respectively.  Evaporation  time  is  longer  than  chemical  reaction  time. 
However,  both  physical  processes  are  necessary  to  achieve  a  higher  degree  of  correlation  for  the 
spark  ignition  model.  While  Figures  7  and  8  present  the  evaporation  and  chemical  reaction 
effects  on  primary  zone  stability  limit,  comparison  of  these  physical  processes  indicates  longer 
chemical  reaction  time  with  little  improvement  in  the  degree  of  correlation  for  the  primary  zone 
flame  stabilization  model. 

The  effects  of  flow  parameter  (Wapi  \/T;j/P.})  and  primary  zone  velocity  on  spark  kernel 
quench  and  volumetric  hot  residence  times  are  shown  in  Figures  9,  10,  11,  and  12,  respectively. 
The  figures  imply  that  spark  ignition  and  primary  zone  flame  stabilization  are  impeded  by 
increases  in  flow  parameter  and  velocity  in  the  primary  zone.  A  significant  finding  is  that 
advanced  engines  with  increased  air  flow  and  reduced  volume  in  the  primary  zone  exhibit  the 
lowest  available  spark  quench  and  flow  residence  time  for  a  specified  flow  parameter.  Therefore, 
they  require  shorter  evaporation  plus  chemical  reaction  time  for  a  successful  engine  start.  The 
adverse  effect  of  an  increase  in  flow  parameter  or  velocity  on  the  initial  phase  of  spark  kernel 
development  (anchored  to  the  electrodes)  can  be  offset  by  an  increase  in  minimum  ignition 
energy.  The  same  adverse  effect  on  primary  zone  flame  stabilization  can  be  offset  by  an  increase 
in  the  primary  zone  volume  and/or  decrease  in  primary  zone  airflow. 

One  of  the  most  significant  parameters  which  influence  both  evaporation  and  chemical 
reaction  time  for  spark  ignition,  and  chemical  reaction  time  for  primary  zone  flame  stabilization 
is  the  root  mean  square  of  the  fluctuating  velocity  component  (u'h  In  this  study,  ur  is  proportional 
to  the  square  root  of  the  combustor  pressure  drop  and  inlet  air  density  (see  Section  II).  The 
beneficial  effect  of  an  increase  in  u'  on  spark  ignition  arises  liecause  evaporation  plus  chemical 
reaction  time  is  reduced  as  shown  in  Figures  13  and  14.  Thus,  less  spark  kernel  quench  time  is 
required  for  successful  ignition.  However,  offsetting  this  advantage  is  the  turbulent  diffusion 
heat  loss  suffered  by  the  spark  kernel  as  it  becomes  detached  and  enters  the  recirculation  zone, 
which  reduces  available  spark  kernel  quench  time.  Figures  15  and  lt>  present  similar  beneficial 
effects  of  u'  and  combustor  pressure  drop  on  evaporation  plus  chemical  reaction  time  for  the 
phenomena  of  flame  stabilization.  The  adverse  effect  of  u'  on  primary  zone  stability  is  due  to  the 
entrainment  of  additional  fresh  material  into  the  recirculation  zone,  which  reduces  available  flow 
residence  time  for  evaporation  and  chemical  reaction. 

An  increase  in  Sauter  mean  diameter  (SMD)  adversely  effects  spark  ignition  and  primary 
zone  flame  stabilization  by  increasing  evaporation  time.  For  higher  values  of  SMD.  the  fuel  is  not 
fully  evaporated.  Therefore,  additional  fuel  is  required  to  ignite  and  maintain  combustion  as 
shown  in  Figure  17. 
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Chemical  Reaction  Plus  Evaporation  Effects  on  Primary  Zone  Stability  Limit  During  Windmill  Ignition 
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Chemical  Reaction  Effects  on  Primary  Zone  Stability  During  Windmill  Ignition 
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Figure  10.  —  Effect  of  Flow  Parameter  on  Volumetric  Hot  Residence  Time  During  Windmill  Ignition 


Figure  11.  —  Effect  of  Primary  Zone  Velocity  on  Spark  Kernal  Quench  During  Windmill  Ignition 
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Figure  13.  —  Effect  of  Root  Mean  Square  Turbulent  Velocity  Effect  on  Windmill  Ignition  Limit 
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Combustor  Pressure  Drop  on  Windmill  Ignition  Limit 


Combustor  Pressure  Drop  Effect  on  Stability  Limit  During  Windmill  Ignition 


Diameter  on  Windmill  Ignition  Limit 


2.  SPOOLDOWN  IGNITION  AND  PRIMARY  ZONE  FLAME  STABILIZATION 


Although  spooldown  and  windmill  are  different  modes  of  ignition,  the  formulation  of  the 
characteristic  time  models  for  spark  ignition  and  primary  zone  flame  stabilization  developed  for 
windmill  ignition  were  also  applicable  to  spooldown  ignition.  In  both  cases,  ignition  must  occur 
at  low  levels  of  combustor  pressure  with  slightly  higher  inlet  temperature  during  spooldown 
ignition.  The  most  significant  difference  between  windmill  and  spooldown  modes  is  the  engine 
operation  prior  to  ignition.  During  windmill  ignition,  constant  airflow  is  provided  to  the 
combustor,  while  during  spooldown  ignition,  the  air  flow  is  continuously  decreasing  until  ignition 
followed  by  engine  acceleration  occurs.  Therefore,  the  spooldown  ignition  data  must  cross  the 
stability  and  ignition  limit  at  the  time  of  ignition  while  the  windmill  ignition  data  can  exist  on 
either  side  of  the  limits. 

The  spooldown  spark  ignition  and  primary  zone  flame  stabilization  models  were  developed 
using  F100(3)  and  PW1128  Development  engine  data  (40  percent  and  25  percent  spooldown). 
Linear  correlation  models  of  spooldown  spark  ignition  and  primary  zone  stability  limits  were 
obtained.  For  the  phenomena  of  spark  ignition,  the  available  spark  kernel  quench  time  has  been 
correlated  against  the  evaporation  plus  chemical  reaction  time  with  a  coefficient  of  determina¬ 
tion  of  0.904  and  standard  deviation  of  3.775  as  shown  in  Figure  18.  Also,  for  the  phenomena  of 
flame  stabilization,  the  available  volumetric  hot  residence  time  has  been  correlated  against 
evaporation  plus  chemical  reaction  time  with  a  coefficient  of  determination  of  0.946  and  standard 
deviation  of  0.187  as  shown  in  Figure  19.  The  models  are  valid  over  the  entire  range  of  conditions 
for  which  data  were  available. 


Separate  plots  illustrating  the  effects  of  evaporation  and  chemical  reaction  time  on  spark 
kernel  quench  time  and  volumetric  hot  residence  time,  are  shown  in  Figures  20  through  23.  Both 
physical  processes  are  necessary  to  achieve  a  higher  degree  of  correlation  for  the  phenomena  of 
ignition  and  stabilization.  However,  evaporation  time  is  longer  for  the  spark  ignition  model  while 
chemical  reaction  time  is  longer  for  the  primary  zone  flame  stabilization  model. 


The  adverse  effects  of  flow  parameter  and  primary  zone  velocity  on  spark  kernel  quench 
time  and  volumetric  hot  residence  time  are  shown  in  Figures  24  through  27.  Similar  results  were 
also  obtained  with  the  windmill  ignition  data.  Future  engines  show  trends  of  decreasing  spark 
ignition  and  stability  limit  with  increasing  flow  parameter  and  primary  zone  velocity  (increasing 
airflow  and/or  decreasing  primary  zone  volume). 


The  adverse  and  beneficial  effects  of  an  increase  in  u'  (increase  pressure  drop  or  inlet 
temperature)  on  spark  ignition  and  flame  stabilization  limits  were  also  observed  with  the 
spooldown  data  as  shown  in  Figures  28  through  31.  Similarly,  for  the  results  obtained  with  the 
windmill  data,  an  increase  in  SMD  requires  higher  equivalence  ratio  for  successful  ignition  as 
shown  in  Figure  32. 
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Chemical  Reaction  Plus  Evaporation  Effects  on  Spooldown  Ignition  Limits 


Evaporation  Effect  on  Spooldown  Ignition  Limit 


Chemical  Reaction  -  msec 
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Parameter  on  Spark  Kernel  Quench  During  Spooldown  Ignition 


Figure  25.  —  Effect  of  Flow  Parameter  on  Volumetric  Hot  Residence  During  Spooldown  Ignition 


ing  Spooldown  Ignition 


Primary  Zone  Velocity  on  Volumetric  Hot  Residence  Time  During  Spooldown  Ignition 


•e  28.  —  Effect  of  Root  Mean  Square  Turbulent  Velocity  on  Spooldown  Ignition 


Combustor  Pressure  Drop  on  Spooldown  Ignition 


Combustor  Pressure  Drop  on  Spooldown  Ignition 


Combustor  Pressure  Drop  on  Primary  Zone  Stability  Limit  During  Spooldown  Ignition 


Bauter  Mean  Diameter  ~  microns 


SECTION  IV 

LEAN  DECELERATION  BLOWOUT  MODEL 


A  lean  deceleration  blowout  model  was  developed  using  F100-PW-200,  F100-ILC  Develop¬ 
ment  and  PW2037  engine  data.  The  blowout  operating  conditions  were  encountered  by 
decreasing  deceleration  fuel  flow  schedules  (i.e.,  decrease  engine  deceleration  time)  and 
consequently  allowing  lower  fuel-air  ratios  and  combustor  pressures  to  occur  (all  limits  were 
deactivated).  The  igniters  were  turned  off  such  that  a  full  combustor  blowout  would  not  relight 
thereby  clouding  the  results.  In  order  to  determine  the  combustor  operating  conditions  during 
blowout,  the  derivative  of  the  combustor  pressure,  with  respect  to  time  was  plotted  as  a  function 
of  the  combustor  pressure  as  shown  in  Figure  33.  As  can  be  seen,  the  blowout  conditions  were 
encountered  at  a  decel  schedule  of  3RU  evidenced  by  the  drastic  change  in  the  derivative  of  the 
combustor  pressure,  while  the  combustor  operated  successfully  at  a  decel  schedule  of  5RU.  The 
corresponding  successful  operating  conditions  were  determined  based  on  an  average  time  interval 
from  the  start  of  the  decel  to  the  time  the  blowout  points  occurred.  The  corresponding  successful 
operating  conditions  were  determined  based  on  an  average  time  interval  from  the  start  of  the 
decel  to  the  time  the  blowout  points  occurred. 


The  essential  feature  of  a  gas  turbine  combustor,  as  far  as  the  lean  deceleration  model  is 
concerned,  is  the  recirculation  zone  which  is  created  in  the  combustor  primary  zone  flow  field 
and  imposed  by  (1)  the  air  entering  through  the  swirl  vanes;  (2)  the  air  jets  from  the  primary 
combustion  holes;  and  (3)  swirler  purge  air,  and  deflector  cooling  air.  One  of  the  important 
functions  of  the  primary  zone  is  to  recirculate  combustion  products  to  mix  and  burn  with  the 
incoming  air  and  fuel.  By  this  means  a  mechanism  of  continuous  ignition  is  established,  and 
combustion  can  be  sustained  over  a  wide  range  of  operating  condit  ions. 

The  key  to  developing  a  successful  model  for  lean  deceleration  blowout  is  the  identification 
of  the  physical  processes  occurring  in  the  flow  field  of  the  primary  zone.  The  model  is  based  on 
the  characteristic  time  approach  and  considers  the  effects  of  three  processes  fundamental  to 
flame  stabilization:  flow  residence,  chemical  reaction  and  droplet  evaporation.  The  region  of 
interest  is  the  shear  layer  between  the  fresh  air  and  the  recirculating  fuel  and  burned  gas 
mixture.  The  turbulent  eddies  present  in  the  shear-layer  region  must  ignite  and  react  before  they 
are  quenched  by  the  relatively  cold  free  stream;  if  they  cannot  ignite,  the  shear-layer  flame  blows 
off.  The  stabilization  models  with  their  corresponding  characteristic  times  presented  in 
Section  II  were  considered  as  candidates  for  the  lean  deceleration  blowout  model.  Comparison  of 
the  degree  of  correlation  for  these  models  indicated  the  selection  of  the  following  characteristic 
times  for  the  lean  deceleration. 

•  For  volumetric  hot  residence  time 


Tvhr 


Vp,P„ 


(msec) 


•  For  evaporation  time 
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•  For  chemical  reaction  time 
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Identification  of  Lean  Decel  Blowout  Points  for  FlOO-ILC  Development  Engine  FX217-21 


An  excellent  correlation  model  of  the  lean  deceleration  blowout  limit  was  obtained  with  the 
aforementioned  engine  data.  As  shown  in  Figures  34  and  35,  the  volumetric  hot  residence  time  in 
the  primary  zone  has  been  correlated  with  the  characteristic  times  of  chemical  reaction  and 
evaporation.  The  relationship  defines  regions  of  stable  and  unstable  combustion.  The  solid  line 
represents  the  lean  deceleration  blowout  limit  over  the  entire  range  of  conditions  for  which  data 
were  available.  The  coefficient  of  determination,  r2  was  0.985. 

The  successful  deceleration  data  points  lie  in  the  region  of  stable  combustion  and  approach 
the  lean  decel  blowout  limit  as  the  deceleration  rate  is  increased  (i.e.,  decreased  engine 
deceleration  time)  as  shown  on  the  dashed  lines  in  Figure  34.  For  a  specified  combustor  geometry 
and  operating  condition  (fixed  volumetric  hot  residence  time),  the  blowout  equivalence  ratio  is 
uniquely  related  to  the  sum  of  chemical  reaction  and  evaporation  time.  Alternately,  an  increase 
in  the  volumetric  hot  residence  time  of  the  primary  zone,  by  reducing  airflow  (decreasing  primary 
zone  velocity)  or  increasing  primary  zone  volume,  will  permit  use  of  lower  equivalence  ratio. 

Plots  illustrating  separately  the  effects  of  evaporation  and  chemical  reaction  time  on  lean 
deceleration  limit  are  shown  in  Figures  36,  37,  and  38.  The  chemical  reaction  rates  are 
significantly  longer  than  the  evaporation  rates  for  well  atomized  fuel  and  low  equivalence  ratio. 
Therefore,  chemical  reaction  time  was  concluded  to  be  the  controlling  process  of  lean 
deceleration  blowout.  Also,  as  shown  by  the  dashed  line  in  Figure  36,  an  increase  in  the 
deceleration  rate  has  no  significant  effect  on  the  evaporation  time,  while  a  decrease  in  the 
volumetric  hot  residence  time  increases  evaporation  time  slightly. 

Figures  39  and  40  present  the  effect  of  the  pressure  drop  and  overall  fuel  air  ratio  on 
chemical  reaction  time,  respectively.  With  increased  pressure  drop  from  one  engine  to  another, 
the  chemical  reaction  time  is  reduced,  since  pressure  drop  induces  mixing  between  hot  gases  and 
fresh  reactants.  However,  for  a  given  engine,  chemical  reaction  time  increases  with  increasing 
deceleration  rate  (i.e.,  decreased  engine  deceleration  time)  until  blowout  occurs  as  shown  in 
Figure  38  The  slight  increase  in  pressure  drop  is  overcome  by  the  decrease  in  flame  speed 
associated  with  the  decrease  in  fuel  air  ratio.  Also,  as  shown  in  Figure  40,  a  decrease  in  combustor 
pressure  results  in  lean  decel  blowout  at  higher  overall  fuel  air  ratio  with  a  longer  engine 
deceleration  time.  Similar  results  were  also  observed  with  the  other  engines. 

The  volumetric  hot  residence  time  is  shown  as  a  function  of  fuel  flow  to  burner  pressure 
ratio  at  blowout  in  Figure  41.  Reduction  in  fuel  flow  to  burner  pressure  ratio  (i.e.,  decreased 
engine  deceleration  time)  increases  volumetric  hot  residence  time.  The  effect  is  less  pronounced 
with  engines  of  lower  flow  residence  time  (increased  front  end  air  flow  or  reduced  primary  zone 
length).  Also,  as  shown  in  Figure  42,  increasing  deceleration  rate  along  a  constant  pressure  line 
with  engines  of  longer  volumetric  hot  residence  time  results  in  a  blowout  at  lower  primary  zone 
equivalence  ratio  and  shorter  engine  deceleration  times. 

The  effects  of  combustor  inlet  temperature  and  pressure  on  lean  deceleration  blowout  limits 
are  shown  in  Figures  43  and  44,  respectively.  Increases  in  inlet  temperature  lead  to  higher 
temperature  combustion  products  and  a  decrease  in  chemical  reaction  time.  An  increase  in 
combustor  pressure  results  in  an  increase  in  flame  speed  and  therefore  a  reduction  in  chemical 
reaction.  Hence,  this  leads  to  a  widening  of  the  stability  limits  with  successful  lean  deceleration 
at  lower  fuel  air  ratio. 
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Figure  36.  —  Evaporation  Effect  on  Lean  Decel  Blowout  Limit 


'igure  39.  —  Effect  of  Pressure  Drop  on  Chemical  Reaction  Time  During  Lean  Decel  Blowout 


'igure  41.  —  Effect  of  Fuel  Flow  to  Burner  Pressure  Ratio  on  Volumetric  Hot  Residence  Time  During  Lean  Decel  Blowout 


Burner  Pressure  Ratio  Versus  Primary  Zone  Equivalence  Ratio 


SECTION  V 


CONCLUSIONS 


Based  upon  the  results  of  the  analyses  conducted  in  this  study,  a  number  of  conclusions 
have  been  drawn  relative  to  the  effect  of  the  important  physical  processes  in  the  combustor  on 
altitude  ignition  and  lean  deceleration  blowout  performance,  and  the  relationship  of  these 
processes  to  combustor  operation  conditions  and  design  variables. 

1.  The  characteristic  time  models  for  spark  ignition  and  flame  stabilization 
were  used  to  correlate  windmill  altitude  ignition  and  lean  deceleration 
blowout  limits  of  current  and  advanced  military  and  commercial  engine 
data.  The  correlations  had  near  zero  intercepts  with  good  statistical 
significance. 

2.  Although  the  evaporation  time  is  longer  than  the  chemical  reaction  time  for 
the  phenomena  of  windmill  altitude  ignition,  consideration  of  both  physical 
processes  provides  a  higher  degree  of  correlation. 

3.  Increased  flow  parameter  adversely  affects  spark  ignition  since  it  reduces 
available  spark  kernel  quench  time.  Advanced  military  engines  with 
increased  air  flow  and  reduced  volume  in  the  primary  zone  for  specified 
flow  parameter  exhibit  the  lowest  available  spark  quench  time. 

4.  Windmill  altitude  ignition  is  impeded  by  increased  airflow  and  reduced 
primary  zone  volume  (i.e.,  increased  velocity  in  the  primary  zone)  for 
current  and  advanced  gas  turbine  combustors.  The  increased  primary  zone 
velocity  reduces  available  spark  kernel  quench  time.  Similar  results  were 
also  observed  with  increased  turbulence  intensity  or  the  fluctuating 
component  of  velocity.  The  spark  quench  time  is  reduced  due  to  the 
turbulent  diffusion  heat  loss  suffered  by  the  spark  kernel  as  it  enters  the 
recirculation  zone. 

5.  The  adverse  effect  of  increased  turbulence  intensity  (i.e.,  decreased 
available  spark  kernel  quench  time)  is  offset  by  the  decreased  evaporation 
plus  chemical  reaction  time  for  successful  ignition  and  acceleration  to 
occur. 

6.  Increased  SMI)  requires  higher  fuel  air  ratio  to  ignite  and  maintain 
combustion. 

7.  The  lean  deceleration  blowout  limit  for  the  engine  data  used  in  this  study  is 
controlled  by  chemical  reaction. 

8.  Longer  engine  deceleration  time  is  required  to  avoid  lean  deceleration 
blowout  with  decreased  primary  zone  volume  or  increased  combustion  zone 
airflow. 

9.  Increased  pressure  drop  has  a  beneficial  effect  on  lean  deceleration  blowout 
limit  by  decreasing  chemical  reaction  time.  However,  offsetting  this 
advantage  is  the  reduction  in  the  primary  zone  flow  residence  time. 
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10.  Decreased  combustor  pressure  and  temperature  adversely  effect  stability 
limits  with  lean  deceleration  blowout  at  higher  fuel  air  ratio. 

11.  The  critical  parameters  of  applying  the  windmill  altitude  ignition  and  lean 
deceleration  blowout  models  are:  (1)  the  empirical  estimates  of  SMD, 
(2)  the  equivalence  ratio  in  the  primary  zone  and  near  the  ignitor  based  on 
fuel  spray  distribution,  and  (3)  the  turbulence  intensity  in  the  combustion 
zone. 
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when  the  mean  flow  velocity  exceeds  the  turbulent  burning  velocity. 
The  weak  extinction  data  are  shown  to  correlate  with  the  same 
parameters  as  conventional  turbulent  burning  velocity  data. 
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